Comparison of Theory and Experiment 
for a One- Atom Laser in a Regime of Strong Coupling 



A. D. Boozer, A. Boca, J. R. Buck, J. McKeever, and H. J. Kimble 

Norman Bridge Laboratory of Physics 12-33, California Institute of Technology, Pasadena, CA 91125 

(Dated: February 1, 2008) 

Our recent paper reports the experimental realization of a one-atom laser in a regime of strong 
coupling Here we provide the supporting theoretical analysis relevant to the operating regime 
of our experiment. By way of a simplified four-state model, we investigate the passage from the 
domain of conventional laser theory into the regime of strong coupling for a single intracavity atom 
pumped by coherent external fields. The four-state model is also employed to exhibit the vacuum- 
. . Rabi splitting and to calculate the optical spectrum. We next extend this model to incorporate 

, the relevant Zeeman hyperfine states as well as a simple description of the pumping processes 

' in the presence of polarization gradients and atomic motion. This extended model is employed to 

' make quantitative comparisons with the measurements of Ref. 01 for the intracavity photon number 

versus pump strength and for the photon statistics as expressed by the intensity correlation function 

^ : 

^ ■ I. INTRODUCTION 

^ : 

' Althou g^h a number of theoretical analyses related to a one-atom laser have appeared in the literature H, H, ^ IE 
. lESIF^ 9, 10. ^tL IT^ H^. Tl*. 15^. IT^ IitI. these prior treatments have not been specific to the parameter range of our 
^ ' recent experiment as reported in Ref. Because of this circumstance, we have carried out theoretical investigations 
in support of our experimental program, and present comparisons of these model calculations with our measurements 
CO in this paper. In Section II we introduce a simplified four-state model that captures the essential features of the 
operation of our one-atom laser in a domain of strong coupling but which avoids the complexity of the full Zeeman 
substructure of the hyperfine levels in atomic Cesium. Sections III and IV then present in turn semiclassical and 
quantum solutions for this four-state model system. By way of a physically motivated transformation for which the 
length of a Fabry-Perot cavity is made progressively shorter, we utilize these results to investigate the continuous 
passage from a domain in which conventional laser theory is applicable into a regime of strong coupling for which the 
_ , full quantum theory is required. We thereby gain some insight into the relationship of our system to prior theoretical 
treatments related to the definition of the laser threshold and to "thresholdless" lasing [l^, [l^ l^fl 123, l^ll ■ The four- 
' state model is further employed to calculate the intracavity photon number versus pump detuning, thereby exhibiting 
^ . the "vacuum-Rabi" splitting for the atom-cavity system 0, US and to compute the optical spectrum of the 
^ ' intracavity field. 

Q^. In Section V we describe the procedures followed to obtain solutions for an expanded model that incorporates 
^ ' the relevant Zeeman substructure for the Cesium atom (32 atomic states), two modes of the cavity with orthogonal 
. J—^ , polarizations, and a simple model to account for the polarization gradients of the optical fields. Comparisons of 
■ the results from quantum jumps simulations based upon this expanded model with our measurements of the mean 
JH intracavity photon number n versus normalized pump intensity x (Figure 3 of Ref. 0) and with our experimental 
determination of the intensity correlation function g'^^^r) (Figure 4 of Ref. [J) are given in Sections V(a) and V(b), 
respectively. 

Our intent here is not to belabor the comparison of our experiment with prior work on niicro-niasers and lasers, 
for which extensive reviews are available [26l '29*. . Instead, our principal goal is to establish quantitative 

correspondence between our measurements and fundamental theoretical models. Having thereby validated the suit- 
ability of the theoretical treatments, we can then use these models to inform further experimental investigations of 
the atom-cavity system. 
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II. FOUR-STATE MODEL 



We begin with a four-state model to describe our experiment in which a single Cesium atom is trapped inside an 
optical cavity as illustrated in Figure ^ Although the actual level structure of the Cesium 65*1/2 ^ 6P3/2 transition 
is more complex due to the Zeeman substructure, this simpler model offers considerable insight into the nature of 
the steady states and dynamics. Following the labelling convention in Fig. ^ we introduce the following set of 




FIG. 1: Illustration of a one-atom laser, (a) The atom is located in a high-Q optical cavity of decay rate k, and is driven by 
the fields f23,4. (b) Inset of the atomic level scheme relevant to our experiment with the 651/2 ^ 6P3/2 transition in atomic 
Cesium. The "lasing" transition is from the excited level F = 3' to the ground level F = 4. Pumping of the excited 3' level 
is by way of coherent excitation from a laser with Rabi frequency Qs. Effective decay from the ground 4 level is provided by 
the combination of a second field with Rabi frequency Q4 and spontaneous decay 4' 3. Various radiative decay rates yij 
appropriate to the D2 line in Cs are given in the text. 



Hamiltonians Hi in a suitably defined interaction picture {h — 1): 

Hi = g43(a^(7g4,e3 + 0-e3,34a), (1) 
H2 = 2^3('5'g3,e3 + O'eS.gs), 

H4 - (AAc + A4)a^a, 

H5 = A3(Te3^e3 + A4(7e4_e4, 

Htot = Hi + H2 + H3 + H4 + Hr-, . 

In a standard convention, the atomic operators are aij = \i){j\ for states with the association of the F ~ 3,4 

ground and the F' — 3', 4' levels with g3,g4, e3, e4, respectively. The Hamiltonian Hi accounts for the coherent 
coupling of the atomic transition e3 <-> gA to the field of a single mode of the cavity with creation and annihilation 
operators (d^,a). The upper state e3 of the lasing transition is pumped by the (coherent-state) field fis, while 
the lower state g4 is depleted by the field fli as described by {H2,H3), respectively. (Hi^H^) account for various 
detunings, including Aac for the offset between the cavity resonance and the e3 <-^- gA atomic transition, A3 for the 
offset between the field ft3 and the g3 ^ e3 transition, and A4 for the offset between the field il4 and the gA ^ e4 
transition. Beyond these interactions, we also account for irreversible processes by assuming that the atom is coupled 
to a continuum of modes other than the privileged cavity mode, and likewise for the coupling of the cavity mode to 
an independent continuum of external modes. 

With these preliminaries, it is then straightforward to derive a master equation for the density operator p for the 
atom-cavity system |3ll l32| in the Born-Markov approximation. For our model system, this equation is 

§ = -i[Htot,p]+j2^- (2) 
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Here, the terms Li account for each of the various decay channels, and are given expUcitly by 

Li = n(2dpa^ — d^ ap — pa^ a), (3) 

L2 — 733(2(Tg3,e3/5(Te3,33 — CTeS.eSP — /5(Te3,e3), 

— 743(2(Tg4,e3PO'e3,34 — ^e3,e3P — P^e3,e3): 

L4 = 734(2(7g3,e4PO'e4,g3 ~ <^e4,e4P ~ e4,e4) , 

L5 = 744(2(7g4,e4P0'e4,g4 — <5'e4,e4/5 — /5(7e4,e4) , 

where the association of each term Li with the decay processes in Fig. ^ should be obvious. Spontaneous decay of 
the various atomic transitions to modes other than the cavity mode proceeds at (amplitude) rate "fij as indicated in 
Fig. while the cavity (field) decay rate is given by k. 

The master equation allows us to derive a set of equations for expectation values of atom ((T,;.j) and field (a) 
operators. One example is for the atomic polarization (fTg4^e3) on the e3 gA transition, namely 



d{ag4. 



e3 



dt 



- [(733 + 743) + « A3] ((Tg4,e3) (4) 

(ri3((Tg4,33) - ^4{o'e4,e3)) 
+1543 i{^e3,e3d) - (lTg4,g4a)) ■ 



A solution to this equation requires not only knowledge of single-operator expectation values (iTij) and (a), but also 
of operator products such as (i5'e3,e3a)- We can develop coupled equations for such products {aijo} but would find 
that their solution requires in turn yet higher order correlations, ultimately leading to an unbounded set of equations. 

Conventional theories of the laser proceed beyond this impasse by one of several ultimately equivalent avenues. 
Within the setting of our current approach, a standard way forward is to factorize operator products in the fashion 

(a, J a) = ((Tij)(a) + ((<T,.ja) - (<T.,.j)(a)) (5) 

with then the additional terms of the form ({aijd) — {aij){d)) treated as Langevin noise. Such approaches rely 
on system-size expansions in terms of the small parameters (l/ng, 1/iVo), where (no. No) are the critical photon 
and atom number introduced in Ref. P| for our one-atom laser. Within the context of conventional laser theory, 
these parameters are described more fully in Ref. |3ll l32j | , while their significance in cavity QED is discussed more 
extensively in Ref. js^. In qualitative terms, conventional theories of the laser in regimes for which {no, No) S> 1 
result in dynamics described by evolution of mean values (o'ij) and (a) (that are of order unity when suitably scaled), 
with then small amounts of quantum noise (that arise from higher order correlations of order (I/tiq: l/-^o) <C 1 ). 

In the following section, wc discuss the so-called semiclassical solutions obtained from the factorization {(Ji.jd) = 
((7ij)(a) neglecting quantmn noise. In Section IV, we then describe the full quantum solution obtained directly from 
the master equation. 

III. SEMICLASSICAL THEORY FOR A FOUR-STATE ATOM 

We will not present the full set of semiclassical equations here since they are derived in a standard fashion from 
the master equation Eq. 121 |32l l3^ . One example is for the atomic polarization {cTg4^e3) on the e3 <-*■ 34 transition, for 
which Eq. 01 becomes 



d{ag4, 



e3) 



dt 



- [(733 + 743) + iAs] ((Tg4,e3) (6) 
-i (03((Tg4,g3) - rJ4(o-e4,e3>) 
-«543 ((0'e3,e3> " g4^g4)) a , 



where a = (a). There is a set of 18 such equations for the real and imaginary components of the various field and 
atomic operators, together with the constraint that the sum of populations over the four atomic states be unity. We 
obtain the steady state solutions to these equations, where for the present purposes, we restrict attention to the case of 
zero detunings I^ac = A3 = A4 = 0. Allowing for nonzero detunings of atom and cavity would add to the complexity 
of the semiclassical analysis because of the requirement for the self-consistent solution for the frequency of emission 
[see, for example, Ref. 35] for the case of a (multi-atom) Raman laser]. 
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FIG. 2; Results from the semiclassical theory as apphed to the atom-cavity system in Fig. 1. (a,b) Intracavity intensity \a\ in 
units of the critical photon number no is plotted as a function of the pump intensity I3 = A threshold for \a\^ is evident 

for 73 ~ 0.8. (c,d) Populations an — (an) versus I3. In (c), population inversion ae3,e3 > <^g4,,g4 occurs over a wide range as the 
pump intensity I3 is increased from 0, including in the threshold region I3 ~ 0.8, with then "population clamping" for CTe3,e3 



as I3 increases beyond threshold. In all cases, the recycling intensity 74 = (j^) — 3 and the detunings Aac = A 



0. 



The semiclassical solutions are obtained for the parameters relevant to our experiment with atomic Cs, namely 

3 17 5 

(733,743,744,734) = (4' 4' Y2' l^-*"^ ' 

where these rates are appropriate to the (amplitude) decay of the levels 6^3/2,^' — 3', 4' — > 65*1/2,7^ = 3,4 with 
'Y — 2tt X 2.6 MHz (i.e., a radiative lifetime r = I/27 — 30.6 ns). The cavity (field) decay rate k is measured to be 
K = 27r X 4.2 MHz. The rate of coherent coupling 543 for the e3 ^ gA transition (i.e., 6P3/2, F' — 3' ^ QS1/2, 7^ = 4) is 
calculated from the known cavity geometry (waist and length) and the decay rate 7, and is found to be 343 = 27r x 16 
MHz based upon the effective dipole moment of the transition. 

Examples of the resulting steady-state solutions for the intracavity intensity |ap together with the populations 
an of the four atomic states are displayed in Figure |2 Parts (a) and (c) of the figure illustrate the behavior of jap 
and au around the semiclassical threshold as functions of the pump intensity 73. Parts (b) and (d) explore these 
dependencies over a wider range in 73. For fixed ratios among the various decay rates as in Eq. EJ the semiclassical 
solutions for |ap/no as well as the various populations an plotted in Fig. |21 depend only on the critical atom number 
A^o (or equivalently, the cooperativity parameter Ci — 1/Nq for a single atom in the cavity). Hence, as emphasized in 
the Supplementary Information published with our paper Ref . |lj , these steady state solutions from the semiclassical 
theory are independent of the cavity length I, and provide a point of reference for understanding "lasing" for a single 
atom in a cavity. This is because Nq = is independent of cavity length I for a cavity with constant mirror 
refiectivity and cavity waist wq. 

Importantly, the semiclassical theory predicts threshold behavior for parameters relevant to our experiment, in- 
cluding inversion cre3,e3 > crg4,g4 in the threshold region, although this is not essential for Raman gain for g3 — > gA 
via e3. One atom in a cavity can exhibit such a "laser" transition for the steady state solutions in the semiclassical 
theory because the cooperativity parameter Ci 3> 1. Indeed, in these calculations we used our experimental value for 
the cooperativity parameter Ci = 1/A^o — 12. Among other relevant features illustrated in Fig. El is the quenching 
of the laser emission around 73 ~ 6.5, presumably due to an Autler-Townes splitting of the excited state e3 at high 
pump intensity 8J. 



5 



A. Relationship to a Raman laser 

In many respects our system is quite similar to a three-level Raman scheme, for which there is an extended literature 
(e.g., Ref. }35 j and references therein). In fact we have carried out an extensive analysis of a Raman scheme analogous 
to our system in Fig. ^ Pumping is still done by the field fis on the 3^3' transition. However, recycling 4 ^ 4' — > 3 
by the field ^4 and decay 734 is replaced by direct decay 4 — > 3 at a fictitious incoherent rate of decay (334 with level 
4' absent. In all essential details, the results from this analysis are in correspondence with those presented from our 
four-level analysis in this section. In particular, the threshold onsets in precisely the same fashion as in Fig. El^a), and 
the output is "extinguished" at high pump levels for fl^. This turn-off appears to be associated with an AC-Stark 
splitting of the excited 3' level by the fl^ field that drives the 3' ^ 4 level out of resonance with the cavity due to the 
splitting of the upper level 3'. Over the range of intensities explored in this section, the "quenching" behavior seems 
to be unrelated to any coherence effect associated with the combination of the field and decay 734. 

IV. QUANTUM THEORY FOR A FOUR-STATE ATOM 

A one-atom laser operated in a regime of strong coupling has characteristics that are profoundly altered from the 
familiar case (described e.g. in Refs. 0,113), for which the semiclassical equations are supplemented with (small) 
quantum noise terms. The question then arises as how to recognize a laser in this new regime of strong coupling, 
where we recall the difficulty that this issue engenders even for systems with critical photon number much greater 
than unity T^, "201 1^ . The perspective that we adopt here is to investigate the continuous transformation of a 
one-atom laser from a domain of weak coupling for which the conventional theory should be approximately valid into 
a regime of strong coupling for which the full quantum theory is required. 

Towards this end, we consider a scenario in which the cavity length (and hence its volume) is gradually reduced 
from a "large" value for which the conventional theory is valid to a "small" value for which the system is well into 
a regime of strong coupling. As illustrated in Figure |3] , this transformation is assumed to be under conditions of 
constant cavity waist wo and mirror reflectivity R, in which case scaling the length by a factor / causes the other 
parameters to scale as follows: 



I - 


- lf = 


fl, 


9 - 


9f = 


9/f'^' 


K - 


Kf ^ 


= 


7 - 


7, 




Nq - 


A^o, 




no - 


-> no/ 


= fno . 



Recall that in the semiclassical theory illustrated in Fig. |21 the quantity jap/no / is invariant under this transformation. 
By contrast, the role of single photons becomes increasingly important as the cavity length is reduced (i.e., uq / becomes 
ever smaller), so that deviations from the familiar semiclassical characteristics should become more important, and 
eventually dominant. 




FIG. 3: Illustration of the scaling transformation considered in Eqn. |H| whereby the length of a spherical mirror Fabry-Perot 
cavity is transformed I fl while the cavity waist wo and the atomic position are held constant. The atom is indicated by the 
"dot" in the center of the cavity mode. 
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FIG. 4: The mean intracavity photon number n/nof (blue) and normalized intensity correlation function ^'^''(O) (red) are 
plotted as functions of pump intensity J3 = (^3/27)^ in (a)-(d). In (a)-(c), the cavity length is made progressively shorter 
(2500^0, lOOio, lo), where = 42.2 /^m is the length of our actual cavity. The corresponding saturation photon numbers are 
no/ =(33.0, 1.32, 0.013). n/no/ and (7'^'(0) are calculated from the quantum theory for the four-state system in Fig. Q while 
\a\'^ /nof given by the black curve is from the semiclassical theory, (d) n (blue), ^'^•'(O) (red), and the Mandel Q parameter 
(green) shown over an extended range of pump intensity J3 for I = Iq. In all cases, I4, — (^4/27)^ = 2, the 3^4' and 4 — > 4' 
transitions are driven on resonance, and the cavity detuning uoca ~ 0. Other parameters are as given in the text. 



A. Field and atom variables for various cavity lengths 

Framed by this perspective, we now present results from the quantum treatment for a four-state model for the 
atom. Our approach is to obtain steady state results for various operator expectation values directly from numerical 
solutions of the master equation given in Eq. |21by way of the Quantum Optics Toolbox written by S. Tan |36|. Since 
such numerical methods are by now familiar tools, we turn directly to results from this investigation presented in 
Figs. EE 

These figures display the behavior of various characteristics of the atom-cavity system as the cavity length is reduced 
from I ~ 25OOZ0 to I = 100^0 to I — Iq to I = ^o/99, where = 42.2/im is the actual length of our cavity. Figure 
0] provides an overview of the evolution and is reproduced from the Supplementary Information in Ref . [ij , while 
Figures I5I9I provide more detailed information about the intracavity field and atomic populations. 

Figure EJa-c) and part (a) in Figs. El IHI and [T] display the mean intracavity photon number n/riQ/ (where no/ 
is calculated for the particular length), and compare this result to |ap/no/ from the semiclassical theory. The 
correspondence is close in Figs. EJa) and [SJa) since riQ/ = 33 in this case, but becomes increasingly divergent in 
Figs. Elb) andEJa) for which np/ = 1.3, and in Figs, ^c) and[3a) for which no/(/ = 1) = ^^o = 0.013 (as in our 
experiment). 

In qualitative terms, the peak in each of the curves for n/riQ/ in Figs. and|7|arises because of a "bottleneck" 
in the cycle g3 —f e3 —^ —^ e4: —^ g3. For our scheme with one atom in a cavity, this cycle can proceed at a rate 
no faster than that set by the decay rate 734. For higher pump intensities I3, the quenching of the emission displayed 
by the semiclassical theory becomes less and less evident with decreasing / as the coherent coupling rate g becomes 
larger in a regime of strong coupling. 

Part (b) in Figs. [S] [HI and [T] shows the populations an of the four states. A noteworthy trend here is the rapid 
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FIG. 5: Steady state solutions as functions of pump intensity Is obtained from the numerical solution of the master equation 
l^for the four-state atom in a cavity illustrated in Fig. Here, the cavity length I = 2500^0, where lo = 42.2/im is the cavity 
length in our experiment, (a) Mean intracavity photon number n normalized to the saturation photon number no/ =33 (in 
blue). The corresponding result for |ap/no/ from the semiclassical theory is given by the black curve, (b) Populations au of the 
four states as labelled, (c) Mean intracavity photon number n (blue), Mandel Q parameter (green), and intensity correlation 
function g^'^\0) (red), (d) Ratio R of photon flux from the cavity mode k/ n as compared to the rate of atomic fluorescence 
743cre3,e3 for the exclted state e3. In all cases, the depleting intensity Ia = (§^)^ = 3 and the detunings Aac = A3 = A4 = 0. 
Field and atom decay rates are as specified in the text. 



reduction of the population (Jes.es, with decreasing cavity length. Again, the rate g becomes larger as I is reduced, and 
eventually overwhelms all other rates, so that population promoted to this state is suppressed. 

Figure 0] and part (c) in Figs. El El and [7| address the question of the photon statistics by plotting the Mandel Q 
parameter (or equivalently the Fano factor F ^ Q -\- 1) as well as the normalized second-order intensity correlation 
function 5(2^(0) [33. As shown in Fig. [^a), for large / = 2500/o^he region around the semiclassical threshold 
displays the familiar behavior associated with a conventional laser jS^, 0, IS IH, 113 , namely that ^'^^•'(O) evolves 
smoothly from ^(^^(O) ~ 2 below the semiclassical threshold to g^'^\Q) « 1 above this threshold. Furthermore, Fig. 
jSJc) shows that the Mandel Q parameter has a maximum in the region of the threshold 19] . Beyond this conventional 
(first) threshold, the Mandel Q parameter in Fig. El (c) also exhibits a second maximum, that has been described as 
a "second" threshold for one-atom lasers Q, and g^'^\0) rises back from 1 to 2. With decreasing cavity length, these 
features are lost as we move into a regime of strong coupling. For example, the two peaks in Q merge into one broad 
minimum with Q < indicating the onset of manifestly quantum or nonclassical character for the emission from the 
atom-cavity system. 

Finally, part (d) in Figs. El El and [7| presents results for the ratio _R, where 

R = i (9) 

7430'e3,e3 

gives the ratio of photon flux nffi from the cavity mode to the photon flux 743cre3,e3 appearing as fluorescence into 
modes other than the cavity mode from the spontaneous decay e3 gA. For a conventional laser, Kffi <^ 7430'e3.e3 
below threshold, and Kffi ^ 743Ce3,e3 above threshold, with the laser threshold serving as the abrupt transition 
between these cases in the manner of a nonequilibrium phase transition |3 [s^ • As illustrated in Fig. {7\ no such 
transition is required in the regime of strong coupling; R ^ 1 from the onset as the pump I3 is increased. This 
behavior is analogous to the "thrcsholdless" lasing discussed in Rcfs. [Tsl 12(1. 121I l2^ and reviewed by Rice and 
Carmichael [l^ . 

For the system illustrated in Figure El the progression in length reduction has a limit at I = Xo/2 corresponding 
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FIG. 6: Steady state solutions as functions of pump intensity /s obtained from the numerical solution of the master equation 
121 for the four-state atom in a cavity illustrated in Fig. Q Here, the cavity length I = lOO/o, where lo = 42.2/im is the cavity 
length in our experiment, (a) Mean intracavity photon number n normalized to the saturation photon number no/ = 1.3 (in 
blue). The corresponding result for \af /nof from the semiclassical theory is given by the black curve, (b) Populations an of the 
four states as labelled, (c) Mean intracavity photon number n (blue), Mandel Q parameter (green), and intensity correlation 
function g^^\0) (red), (d) Ratio R of photon flux from the cavity mode k/ n as compared to the rate of atomic fluorescence 
7430"e3,e3 for the excited state e3. In all cases, the depleting intensity I4 = (^)^ = 3 and the detunings Aac = A3 = A4 = 0. 
Field and atom decay rates are as specified in the text. 

to a Fabry-Perot cavity with length equal to the lowest order longitudinal mode Xq/2, where Aq = 852.3 nm is the 
wavelength of the cavity QED transition. To reach this limit from the length Iq appropriate to our actual cavity, we 
must scale Iq flo with / = 1/99. In a continuation of the sequence shown in Figs. 0,13 and[7| we display in Fig. 
results for such a cavity with I — Ao/2. Note that although Ci = 1/Nq ~ 12 is invariant with respect to this scaling and 
the saturation photon number is reduced to ng/ = 1.31 x 10~^, nevertheless the atom-cavity system has passed out of 
the domain of strong couphng, even though (no/, A'^o) <C 1. This is because strong coupling requires that go ^ (7, k), 
so that (no,7Vo) ^ 1 is a necessary but not sufficient condition for achieving strong coupling. For the progression 
that we are considering with diminishing length (but otherwise with the parameters of our system), I = Ao/2 does 
not lie within the regime of strong coupling (343/7 = 61,343/^ = 0.40), but rather more toward the domain of a 
"one-dimensional atom", for which k 3> g^/^ ^ 7 (see, for example, Refs. |49,l5i| for theoretical discussions and a 
previous experimental investigation). In this domain of the Purcell effect 26^ 28, 29, 30], the fractional emission into 
the cavity mode as compared to fluorescent emission into free space for the 3' — > 4 transition is characterized by the 
parameter 

2(^^(43) 

/343 = ^ ^ 0.99, (10) 

where c|^^^ = Ci x (7/743) 48. 

As compared to Figs. El (HI and [71 a noteworthy feature of the regime depicted in Fig. |S1 is the absence of a 
dependence of ^''^-'(O) on the pump level 13. In fact, g^^-'(O) ~ over the entire range shown, so that the cavity field is 
effectively occupied only by photon numbers and 1. In correspondence to this situation, the Mandel Q parameter in 
Fig. IHIc) is essentially given by the mean of the intracavity photon number, Q ~ — n, with fi <C 1. Furthermore, the 
dominance of emission into the cavity mode over fluorescence decay becomes even more pronounced than in Fig. EJd) , 
as documented by the ratio R in Fig. |Htd). In agreement with ex pec tation set by Eq. ^1 note that R ~ f343 / {1 — . 
All in all, the "bad-cavity" limit specified by k ^ g^/i^ ^ 7 (toward which Fig. |Slis pressing) is a domain of 
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FIG. 7: Steady state solutions as functions of pump intensity Iz obtained from the numerical solution of the master equation 
121 for the four-state atom in a cavity illustrated in Fig. Q Here, the cavity length I = lo, where lo ~ 42.2/^m is the cavity 
length in our experiment, (a) Mean intracavity photon number n normalized to the saturation photon number no = 0.013 (in 
blue). The corresponding result for |ap/no from the semiclassical theory is given by the black curve, (b) Populations an of the 
four states as labelled, (c) Mean intracavity photon number n (blue), Mandel Q parameter (green), and intensity correlation 
function ^'•^•'(O) (red), (d) Ratio R of photon flux from the cavity mode k n as compared to the rate of atomic fluorescence 
7430"e3,e3 for the excited state e3. In all cases, the depleting intensity I4 = (^)^ = 3 and the detunings Aac = A3 = A4 = 0. 
Field and atom decay rates are as specified in the text. 



single-photon generation for the atom-cavity system, which for / <C 1 has passed out of the regime of strong coupling. 

Figures 13 El [7| and |S1 provide a step-by-step description of the evolution of the atom-cavity system from the domain 
of conventional laser theory {I ^ Iq as in Fig. [31 with / = 2500), into the regime of strong coupling {I = Iq as in Fig. 
Ulwith / = 1), and then out of the strong-coupling regime into the Purcell domain {I = /o/99 ~ Ao/2 as approached 
in Fig. [HI with / = 0.01) [2^ [s^ • We now attempt to give a more global perspective of the scaling behavior of 

the atom-cavity system by examining various field and atomic variables directly as functions of the scale parameter 
/ = l/lo- A particular set of such results is displayed in Figure where the pump intensity = 3 is fixed near the 
peak in the output from the semiclassical theory in Fig. [21 and the recycling intensity is held constant at = 3. 

In Fig. inja) the mean intracavity photon number n is seen to undergo a precipitous drop as the cavity length 
is made progressively shorter (i.e., increasing /~^^^, since I oc /). However, when n is normalized to the critical 
photon number ng/ i the quantity n/nof is seen to approach unity for small /^^^^ (i-e., long cavities with I 3> lo) as 
appropriate to the conventional theory in Fig. jSJ. With increases in /~^/^ (i.e., shorter cavity lengths), n/riQf rises to 
a maximum around / ~ 3 for strong coupling with I Iq as in Fig. [71 before then decreasing to approach a constant 
value for yet larger values of f~^l'^ as the system exits from the domain of strong coupling. 

Also shown in Fig. [Hl^a) are the quantities g^^-' (0) and Q + \ that characterize the photon statistics of the intracavity 
field. As previously noted, (?^^''(0) lies in the range 1 < g^^^O) — 2 for conventional laser theory, but drops below 
unity in the regime of strong coupling and approaches zero for / <C 1. In this same limit of very small cavities in the 
Purcell regime, Q ~ — n. 

Fig. [HTb) displays the populations for the four-state system as functions of f^^^^. For the conventional regime 
with /^^'^ 1 , there is population inversion, cre3,e3 > <^gi.g4 (which was shown in Fig. [21for small values of but 
this possibility is lost for increasing f~^^^ (i.e., decreasing cavity length). Strong coupling dictates that the rate g 
dominates all others, so that appreciable population cannot be maintained in the state e3. Finally, Fig. [Hid) displays 
the dependence of the ratio R = (K/n)/(743cre3_e3) on f^^^^. From values i? < 1 in the conventional domain, R rises 
monotonically with decreasing cavity length reaching the plateau i? 3> 1 specified by Eq. ^1 
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FIG. 8: Steady state solutions as functions of pump intensity Is obtained from the numerical solution of the meister equation 
|5|for the four-state atom in a cavity illustrated in Fig. Here, the cavity length / = Zo/99 — Ao/2 (i.e., / = 1/99), where 
lo = 42.2/im is the cavity length in our experiment and Ao ~ 852.3 nm is the wavelength of the cavity QED transition, (a) Mean 
intracavity photon number n normalized to the saturation photon number no/ = 1-31 x 10"'* (in blue). The corresponding 
result for |a|'^/no/ from the semiclassical theory is given by the black curve, (b) Populations an of the four states as labelled. 

(c) Mean intracavity photon number n (blue), Mandel Q parameter (green), and intensity correlation function g^^\0) (red). 

(d) Ratio R of photon flux from the cavity mode k/ n as compared to the rate of atomic fluorescence "1430-^3, es for the excited 
state e3. In all cases, the depleting intensity I4, = (§^)^ = 3 and the detunings Aac = A3 = A4 = 0. Field and atom decay 
rates are as specified in the text. 



B. Vacuum-Rabi splitting 



In the preceding discussion, we have compared various aspects of our one-atom system with conventional lasers and 
have restricted the analysis to the case of resonant excitation with A3 = 0. Our actual system operates in a regime 
of strong coupling, so that there should be an explicit manifestation of the "vacuum-Rabi" splitting associated with 
one quantum of excitation in the 4^3' manifold [2^ |3 @| • 

To investigate this question, we consider the dependence of the average intracavity photon number n on the detuning 
A3 of the pump field ^Is, with the result of this analysis illustrated in Fig. ^| For weak excitation I3 ^ 1 (well 
below the peak in Fig. |Z|^a)), the intracavity photon n is maximized around A3 — ±343 (and not at A3 — 0) in 
correspondence to the eigenvalue structure for the gA e3 manifold in presence of strong coupling. The excited state 
eS is now represented by a superposition of the nondegenerate states \ip±) whose energies are split by the coupling 
energy ±%43 . However, for large pump intensities I3 ~ 10, this splitting is lost as the Autler-Townes effect associated 
with the pump field on the g3 <-> e3 transition grows to exceed 1743. 



C. Optical spectrum of the cavity emission 

A central feature of a conventional laser is the optical spectrum of the emitted field, defined by 

/+CX; 
dT{\im{a\t)a{t + T))}exp{-inT) , (11) 
-00 

where as in Eq. n (at, a) are the creation and annihilation operators for the single-mode field of the cavity coupled 
to the atomic transition e3 gA. The results for the Schawlow-Townes linewidth are well-known and will not be 
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FIG. 9: Scaling behavior of various quantities as the cavity length I is varied, where / — l/lo, and Iq — 42.2 for our actual 
cavity. Note that g cc f^^^^ and k cc /^^ and that the range in / corresponds to that spanned by Figs. |3to|H| namely 
0-01 ^ / 2500. (a) Mean intracavity photon number n (blue), the Mandel Q parameter (Q + 1 in green), and the intensity 
correlation function g^^\0) (red), (b) Populations an of the four states as labelled, (c) Mean intracavity photon number n 
normalized to the saturation photon number no/ — nof — 0.013 x /. (d) Ratio R of photon flux from the cavity mode K/n as 
compared to the rate of atomic fluorescence 743cre3,e3 for the excited state e3, where «/ = k/ f. In all cases, the pumping and 
recycling intensities 13,4 = 3 and the detunings Aac = A3 — A4 — 0. Field and atom decay rates are as specified in the text. 



discussed here [sil Is^. Is^. IssL Is^ . Instead, in Fig. ^]we present results specific to the domain of operation of our 
system. 

For the choice of parameters corresponding to Fig. [7| ^{^) in Fig. Illf a) exhibits a pronounced two-peak structure, 
with the positions of the peaks corresponding to the Autler-Townes splitting of the ground state by the recycling field 
n^. Contrary to what might have been expected from the analysis of the previous section, $(0) shows no distinctive 
features associated with the vacuum-Rabi splitting of the excited state. For reduced values of pumping and recycling 
intensities 13^4 = 0.5, there are small features in the optical spectrum at f2 « ±543, as is illustrated in Fig. 1111 when 
is plotted on logarithmic scale. With respect to the complex degree of coherence [s^l, the coherence properties 
of the light from the one-atom laser in the regime of strong coupling are set simply by the inverse of the spectral 
width of which can be determined from the plots in Fig. ITTI 

The curves shown in Fig. 1111 are calculated by way of the quantum regression theorem applied to the four-state 
system of Fig. ^ From the quantum regression theorem, we have that the two time correlation function in Eq. Illl is 
given by 

(at(O)a(T)) = Tr[pssa^{0)a{T)] = Tr[p{T)a{Q)] 
where p(t) is obtained by numerically evolving 

Pa = PssaHO) 

under the master equation, and pss is the steady state density matrix. By Fourier transforming the correlation 
function according to Eq. 1111 we obtain the optical spectrum. 

The optical spectrum of the emitted light from our cavity could in principle be measured by way of heterodyne 
detection. The cavity output would be combined on a highly transmissive beam splitter with a local oscillator beam 
that is frequency shifted by an interval Alo that is large compared to the range of frequencies in the output field. 
The optical spectrum is then obtained by taking the Fourier transform of autocorrelation function of the resulting 
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FIG. 10: The mean intracavity photon number n versus the detuning A3 (in cycles/sec) of the pump field ^3, where A3 = 
corresponds to the transition frequency 0^33. The three curves are for increasing pump intensity (i) I3 = 0.1, (ii) I3 — 1.0, (iii) 
I3 = 10.0. The arrows indicate the positions of the expected "vacuum- Rabi" peaks at ±343, where gi43/27r — 16 MHz. In all 
cases, the recycling field Q4, is on resonance A4 = and has intensity I4 — 3. 
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FIG. 11: The optical spectrum ^{v) as a function of frequency offset 1/ (in cycles/sec, Cl = 2m/), where v = Q corresponds to 
the transition frequency uja-3- (a) Three spectra '^{v) for increasing pump intensity (i) I3 = 0.1, (ii) I3 = 1.0, (iii) I3 = 10.0, 
with the recycling intensity 74 = 3 in all cases. The overall normalization of ^{v) is arbitrary, but is common for the three 
cases, (b) 4>(i^) on a logarithmic scale for decreased intensities I3 — I4, — 0.5, with the peak value of $ scaled to unity. The 
arrows indicate the position of the expected "vacuum- Rabi" peaks at ±1743, where (j43/27r = 16 MHz. In all cases in (a) and 
(b), the pumping field fl^ and the recycling field fl^ are on resonance with their respective transitions (A3 = = A4). 

heterodyne current. Although we have not carried out this procedure experimentally, it is straightforward to model 
using a quantum jumps simulation of the four state model. We have computed such spectra for several values of 
I3, using a local oscillator flux equal to k. This is an experimentally reasonable value, since it is small enough so 
as to not saturate the detectors, yet large enough that, as our further simulations indicate, increasing the flux does 
not significantly change the resulting spectrum. The results for the spectrum obtained from this quantum jumps 
simulation agree reasonably well with results from the quantum regression theorem presented in Fig. 1111 



13 



V. QUANTUM THEORY INCLUDING ZEEMAN STATES AND TWO CAVITY MODES 

In an attempt to provide a more detailed quantitative treatment of our experiment, we have developed a model 
that includes all of the Zeeman states {F,mp) for the F = 3,4 ground levels and the F' = 3', 4' excited levels of 
the 6Si/2 6P3/2 transition in atomic Cesium, of which there are 32 in total. We also include two cavity modes 
with orthogonal linear polarizations to describe the two nearly degenerate TEMqo modes of our cavity |43|, with 
three Fock states for each mode {|0), |1), |2)}. The total dimension of the Hilbert space for this set of atomic and 
field states is then d= 32 x3x3 = 288, making it impractical to obtain steady state solutions from the master 
equation directly. Instead, we employ the Quantum Optics Toolbox js^ to implement a quantum jumps simulation, 
with various expectation values computed from the stochastic trials. 

In broad outline, our expanded model includes Hamiltonian terms of the form of Eq. ^ with now the terms atj 
generalized to incorporate each of the various Zeeman states. Likewise, the coherent coupling of the atom to the 
cavity takes into account two orthogonally polarized modes (a, 6). The operators Li are similarly modified to obtain 
a new master equation that includes the full set of decay paths among the various states (i.e., cr±, tt transitions), as 
well as the associated quantum collapse terms in the simulation. 

We attempt to describe the dynamics arising from the complex state of spatially varying polarization associated with 
the $^3^4 beams by way of the following simple model. In a coordinate system with the x, z directions perpendicular to 
the cavity axis along y, the 03^4 beams propagate along x, z with orthogonal o± configurations. The helical patterns 
of linear polarization from pairs of counter-propagating beams then give rise to terms in the interaction Hamiltonians 
-^2,3 of the form 

^2 = ^r!3[(S^3,e3 + S^3,93)sin(^3.) (12) 



-(S^3,e3 + Sf3,s3)sin(^3.)] 

1 
"2 



-fi3[(S^3,e3 + S^3,S3)(COS(03X) + COs(03.))] 



and similarly for iJ3 to describe the f24 beams with independent phases (Q/^x, diz)- Here fis and are Rabi frequencies 
corresponding to the incoherent sum of the intensities of the four individual beams. In Eq. 1121 the operators S^g^gj 
are linear combinations of various atomic projection operators for the diverse Zeeman-specific transitions for linear 
polarization along x, y, z, and are given explicitly by 

^g3,e3 ~ ^ "y^ (^gs'^.eS ~ ^g3'",e3) (13) 



^ 

^g3,e3 = ^s3,e3 (1^) 



where 



^3,e3 - 2^ 2^ |53, m) (3, m; 1, g|4, m') {gA, m'\ (16) 



The phases 9i arise from the spatial variations of the polarization state of the Vl^^a beams, and are given, for example, 
by O^x = k^xX with k^x as the wave vector of the pair of 0.^ beams propagating along x. 

The ^3^4 beams tend to optically pump the atom into dark states, with this pumping counterbalanced by atomic 
motion leading to cooling 44] and by any residual magnetic field. In our case, imperfections in the FORT polarization 
[4^ [4^ result in a small pseudo-magnetic field along the cavity axis y 45] with peak magnitude By ~ 0.75 G. This 
pseudo-field By is included in our simulations and tends to counteract optical pumping by the VI^a beams into dark 
states for linear polarization in the x — z plane, 9^x = ^3z ^4.t = O^z = 7r/2, but has no effect for polarization along 
the cavity axis y, Osx = Osz = d^x = 6*4^ = 0. 

Overall, the operation of our driven atom-cavity system involves an interplay of cycling through the levels g3 — > 
e3 — > (74 — > e4 — > g3 to achieve output light on the e3 — *■ gA transition, and of polarization gradient cooling for 
extended trapping times. This latter process involves atomic motion through the polarization gradients of the $^3,4 
beams and is greatly complicated by the presence of By . The detunings and intensities of the ^^3^4 beams are chosen 
operationally such as to optimize the output from our one-atom laser in a regime of strong coupling, while at the 
same time maintaining acceptable trapping times, as shown in Fig. 2 of Ref. jl,]. 
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A. Mean intracavity photon number as a function of pump intensity 

In this section, we present simulation results for the mean intracavity photon number versus pump intensity. In 
qualitative terms, we should expect that the output flux ku predicted from the full multi-state model is significantly 
below that calculated from the four-state model presented in Section IV. This is because the atom necessarily spends 
increased time in manifolds of dark states associated with the pumping by the fls^i beams. 

We can modify the four level model to account for these effects by reducing the decay rate 734 734. The slower 
cycling of the atom due to the reduction of 734 approximates, in a phenomenological way, the slowing effect on the 
recycling of the atom due to optical pumping into dark states. We find that a value 734 = 0.07 x 734 gives a good fit 
to the data (Fig. I12r a.b)). We plot the intracavity photon number versus x = {7 14), since we estimate that 
either measured intensity alone is uncertain by a factor of about 2, but the ratio is known much more accurately. 

For the multi-level simulation, we use two different models to generate mean intracavity photon number versus 
pump intensity curves. In the first model, we neglect the motion of the atom and attempt to capture the essential 
features of the optical pumping processes via a single constant phase = O^x = dsz = Oax = ^4z ■ The choice = Q gives 
no output light, since the f23^4 beams pump the atom into dark states. The value 9 = it 12 chosen for the comparison 
in Fig. I12r c.dl gives good correspondence between the simulations and our measurements with the adjustment of no 
other parameters. For this curve, we plot the average [fia + fib)/ 2 of the intracavity field for the two cavity modes a 
and h. 

As a second, more sophisticated model, we assume that the atom moves at a constant velocity in the radial direction. 
This gives time dependent phases; for example, if we assume that the x coordinate of the atom is 

x{t) — Xq + Vxt 

then 

^■ix{t) — k^xX — 93x,0 + I^Sxt 

where 03x,o = k^xXo, uj^x = k^xVx- For a single simulation run we randomly choose the velocity of the atom and initial 
phases of the ^z^i pumping beams; the intensities from 20 such runs are averaged for each value of x. The velocities 
are chosen uniformly in the range 10 — 20cm/s, which gives angular frequencies in the range 27rl00 — 200kHz. The 
resulting input/output curve is plotted in Fig ll2r e.f'). As before, we plot the average of the intracavity field for the 
two cavity modes. 

We make no claim for detailed quantitative agreement between theory and experiment, as the simulations are 
sensitive to the parameters which are known only approximately, such as the intensity of the ^3^4 pumping beams 
and the magnitude of the pseudo and real magnetic fields. Also, the simulations neglect a number of features of the 
real system, such as atomic motion in the axial direction, the dependence of the cavity coupling g on the position of 
the atom, and a possible intensity imbalance in the ^Is^i pumping beams. However, the simulations do support the 
conclusion that the range of coupling values g that contribute to our results is restricted roughly to O.bgo ^ 9 ^ go- 
Furthermore, the simulations yield information about the atomic populations, from which we deduce that the rate of 
emission from the cavity nn exceeds that by way of fluorescent decay 3' 4, 743/ ((T3/3'), by roughly tenfold over the 
range of pump intensity I3 shown in Fig. 112b . 



B. Photon statistics as expressed by the intensity correlation function g^^''{r) 

In addition to measurements of n versus pumping rate, we have also investigated the photon statistics of the light 
emitted by the TEMqo mode of the cavity by way of the two single-photon detectors D12 illustrated in Fig. 1 of Ref. 

. From the cross-correlation of the resulting binned photon arrival times and the mean counting rates of the signals 
and the backgroundjwe construct the normalized intensity correlation function (see the Supplementary Information 
accompanying Ref. 

^ ' (: I{t) :)2 ' ^ ' 

where the colons denote normal and time ordering for the intensity operators 7 [STj. 

Two measurements for (7^^-'(t) from Figure 4 of Ref. P| are reproduced in (a,b) of Figs. 1131 and lT^ together with 
results from our quantum jumps simulation from the constant phase model with 9 = 7r/2, in (c,d). In Fig. 1131 we 
again have I4 ~ 13 and the pump intensity I3 is set for operation with x ~ 0.17 near the "knee" in n versus x, 
while in Fig. 1141 the pump level is increased to x ~ 0.83. These measurements demonstrate that the light from the 
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FIG. 12: Comparison of theory and experiment for the intracavity photon number n as a function of pump intensity x = 
{7/9){l3/ 14) for fixed I4 = 13 (corresponding to a measured intensity of 50mW/crn? ). The measurements (points with error 
bars) are from Figure 3 of Ref. (a,b) n versus pump strength x for the four level model with 743 = 0.07 x 743. (c,d) n versus 
pump strength x for the constant phase model with 9 = 7r/2. (e,f) n versus pump strength x for the constant velocity model 
described in the text. The immediate onset of emission supports the conclusion of "thresholdless" lasing. Two independent sets 
of measurements are shown as the red and blue points, and agree reasonably well with each other. Details of the measurements 
can be found in Ref. 1], while the parameters for the simulation are given in the text. 



atom-cavity system is manifestly quantum (i.e^nonclassical) and exhibits photon antibunching g^"^^ (0) < g^^^ (r) and 
sub-Poissonian photon statistics g^^\0) < 1 [33. In agreement with the trend predicted by the four-state model in 
Fig. EJc) (as well as by the full quantum jumps simulation), g^'^\0) increases with increasing pump intensity, with 
a concomitant decrease in these nonclassical effects. The bottleneck associated with the recycling process leads to 
this nonclassical character, since detection of a second photon given the first detection event requires that the atom 
be recycled from the F = A ground state back to the F — 3 ground. In this regard, we point to the prior work on 
pump-noise suppressed lasers in multi-level atomic systems, as for example, in Ref. [471. 

In more quantitative terms, theoretical results for g'-^-' (t) from the full quantum jumps simulation are given in parts 
(c,d) of Figs. 1131 and ITll for x — 0.17 and x — 0.83. The excess fluctuations 5*-^-'(t) > 1 extending over t ~ ±1 /is 
appear to be related to the interplay of atomic motion and optical pumping into dark states 44], as well as Larmor 
precession that arises from residual ellipticity in polarization of the intracavity FORT ,4^ .45j . 
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These results for g'^^ (r) provide a perspective on the issue of whether the cavity is effectively "empty" since n is 
quite small. Based upon the mean photon flux from the cavity, this is a reasonably inference, but it is also misleading. 
The nonzero values for g^^^'ir = 0) ~ 0.3, 0.6 ^ 0.01 in Figs. 1131 and IT^ are in fact due to the presence of more than 
one photon in the cavity. Although the mean intracavity photon number is only h ~ 0.005, this number is comparable 
to the saturation photon number uq ~ 0.013. Indeed, the quantum statistical character of the intracavity field is 
determined from the self-consistent interplay of atom and cavity field as in standard laser theories, even though it 
might appear as this interplay is not relevant to the determination of a dynamic steady state. Figure |5| attempts to 
illustrate this point by investigating the passage from the domain of conventional laser theory through the regime of 
strong coupling and thence into a domain of single photon generation with g^'^^ir = 0) ~ over the entire range of 
pumping conditions. 




FIG. 13: The intensity correlation function g'^' (r) of the one-atom laser. (a,b) g^'^^ (r) for a; ~ 0.17 as experimentally determined 
in Ref. 0|. (c,d) Theoretical result for (?'^-'(t) for x = 0.17 from a quantum jumps simulation with 6 = 7r/2. All traces have 
been "smoothed" by convolution with a Gaussian function of width cr = 5 ns. 



C. Discussion of possible coherence effects 

In Section IIIA we briefly described our analysis of an equivalent Raman scheme to address the question of possible 
coherence effects associated with the Vt^ recycling beam. Beyond this analysis, we have also considered the possibility 
that various other coherent processes associated with the pump fields might be important. One concern relates to the 
possibility that 4- wave mixing processes could be important, as for example, in a wave-mixing process that cycles the 
atom 3— s-3'— >4^4'— >3 46]. From an operational perspective, if there were to be a correlated process involved in 
the cycling of the atom 3-^ 3' ^4-^ 4' ^3, then two photons would be emitted into the cavity mode (the "signal" 
on the 3' ^ 4 transition and the "idler" on the 4' ^ 3 transition). In this case since we employ no filter to block 
the "idler" field separated by 9.2GHz, the measured intensity correlation function g^^^r) for the emitted light from 
the cavity would exhibit bunching around t — 0, instead of the observed antibunching and sub-Poissonian character. 
The measured character of g^'^^r) therefore argues against a coherent process that cycles the atom from an initial 
quantum state and back to that state by way of coherent processes involving coupling to the cavity field. 
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FIG. 14: The intensity correlation function g'^' (r) of the one-atom laser. (a,b) g^^^ (r) for x ~ 0.83 as experimentally determined 
in Ref. Q. (c,d) Theoretical result for g'^'^^r) for x = 0.83 from a quantum jumps simulation with 9 = -n 12. All traces have 
been "smoothed" by convolution with a Gaussian function of width cr = 5 ns. 

We also note that the coherent coupling of the cavity field and atom for the 4' ^ 3 transition is greatly suppressed 
due to the large detuning ~ 9.2GHz, leading to an effective coupling coefficient (ge/z/STr ~ 30 kHz) ^ (7||/27r ~ 5.2 
MHz). Therefore, for whatever mixing processes, the coupling to the external vacuum modes characterized by the 
rate 7|| should dominate that due to ge//- In this regard, note that we have included the effect of off-resonant coupling 
of the 4' excited state in our simulations (which is only ~ 200 MHz detuned) . The relevant process is then excitation 
4^4' via the 1^4 pump field, followed by emission into the cavity mode due to the coherent coupling of the transition 
4' 4. This coupling increases the intracavity photon number by only about 10%, suggesting that coupling for the 
4' ^ 3 transition 9.2 GHz away is negligible. 

In support of these comments, our detailed numerical simulations agree sensibly well with the observed behavior of 
g^'^^ij) (as in Figures El and , and do not include any "wave-mixing" effects. This statement is likewise valid for 
the dependence of photon number versus pump level r2§. Furthermore, as previously discussed, the model calculation 
for a four-state system agrees well in its essential characteristics with a three-state system where the decay of the 
ground state 4 ^ 3 is via an ad hoc spontaneous process (as in a Raman laser) rather than by pumping 4^4' and 
decay 4' — > 3. 

A final general comment relates to the nature of phase-matching (e.g., as applied to 4- wave mixing and parametric 
down conversion) for a single atom in a cavity. For a sample of atoms (or a crystal), there is a geometry that defines 
directions for which fields from successive atoms might add constructively for various waves (e.g., pump, signal, 
idler). Cavities can then be placed around these directions to enhance the processes (e.g., the threshold for an optical 
parametric oscillator is reduced by a factor of the square of the cavity finesse for resonant enhancement of both signal 
and idler fields). Clearly a cavity would be ineffective if its geometry did not match the preferred geometry defined by 
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the sample and pump beams. However, for a single atom as in our experiment, these considerations do not apply in 
nearly the same fashion. The relevant issues are the coherent coupling coefficients gij of the various atomic transitions 
to the cavity field. 

VI. SUMMARY 

We have presented a simplified four- level model which describes the qualitative features of our experiment. We 
have shown how decreasing the cavity length causes the model system to move from a regime of weak coupling, where 
the semiclassical laser theory applies, into a regime of strong coupling, where quantum deviations become important. 
The four-state model predicts many of the observed features of our experimental system, including the qualitative 
shape of the intracavity photon number versus pumping intensity curve, and photon antibunching. 

In addition, to predict quantitative values for comparison with our experimental results, we have developed a 
full multi-level model which correctly describes optical pumping and Larmor precession effects within the Zeeman 
substructure. We have shown that these effects play an important role in describing the observed input/ output 
characteristics of the system, and that by including a simple model for the motion of the atom we can obtain 
reasonable agreement with the experimentally observed curve. We have also used the simulation to calculate intensity 
correlation functions, and have compared these results to measurements of ^^^^(t) from our experiment. 
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